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Abstract 

With the help of a simple variational procedure it is possible to convert the partial sums of order of many divergent 
series expansions f{g) - Yj°^=q '^ng" into partial sums Y!^^={) bng '^", where < w < 1 is a parameter that parametrizes 
the approach to the large-g limit. The latter are partial sums of a strong-coupling expansion of fig) which con- 
verge against f{g) for g outside a certain divergence radius. The error decreases exponentially fast for large A^, like 
g-const.xw " present a review of the method and various applications. 
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1. Introduction 

Variational techniques have a long history in theoretical physics. On the one hand, they serve to find equations of 
motion from the extrema of actions. On the other hand they help finding approximate solutions of physical problems 
by extremizing energies. In quantum mechanics, the Rayleigh-Ritz variational principle according to which the ground 
state energy of a system is bounded above by the inequality 



Eo< J d^xiff*(x)Hi//{x) (1) 

has yielded many useful results. In many-body physics, the Hartree-Fock method has helped understanding electrons 
in metals and nuclear matter. In quantum field theory the effective action approach ^ has contributed greatly to the 
theory of phase transitions. In particular the higher effective actions pioneered by Dominicis [2]. 

A variational method was very useful in solving functional integrals of complicated quantum statistical systems, 
for instance the polaron problem Here another inequality plays an important role, the Jensen-Peierls inequality, 
according to which the expectation value of an exponential of a functional of a functional is at least as large as the 
exponential of the expectation value itself: 

(e-^) > e-<^>. (2) 

This technique was extended in 1986 to find approximate solutions for the functional integrals of many other 
quantum mechanical systems 

An important progress was reached in 1993 by finding a way of applying the technique to arbitrarily high order 
list]. The technique was developed furher in the textbook (o]. This made it possible to perform the approximate calcu- 
lation to any desired degree of accuracy. In contrast to the higher effective action approach, the treatment converged 
exponentially fast also in the strong-coupling limit [TJ. 
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Figure 1: Experimental data of space shuttle experiment by Lipa et al. 



The zero-temperature version of this technique led to a new solution of an old problem in mathematical physics, 
that the results of many calculations can be given only in the form of divergent weak-coupling expansions. For 
instance, the energy eigenvalues £ of a Schrodinger equation of a point particle of mass m 



moving in a three-dimensional potential 



can be given as a series ing/ar 



q2 

— + y(x) 

2M 



V(x) 



i//(x) - Ei//(x) 



-X + gx 



n=0 



(3) 



(4) 



(5) 



The coefficients a„ grow exponentially fast with n. The series has a zero radius of convergence. For the ground state 
it reads 

There exist similar divergent expansions for critical exponents which may be calculated from weak-coupling 
expansions of quantum field theories and are experimentally measurable near second-order phase transitions. One 
of these is the exponents a which determines the behavior of the specific heat of superfluid helium near the phase 
transition to the normal fluid. It has been measured with extreme accuracy in a recent satellite experiment [8J. The 
result agrees very well with the value of the series for a as a power series in g/m in the strong-coupling limit m ^ 

iH 

In many more physical examples the properties are found by evaluating divergent weak-coupling series in the 
strong coupling limit. 

In this lecture I shall present the main ideas and sketch a few applications of Variational Perturbation Theory. 



2. Quantum Mechanical Example 

In order to illustrate the method let us obtain the strong-coupling value of the ground state energy We introduce 
a dummy variational parameter by the substitution 

oj^ + (co^ - QZ) = ^Q2 + gr, (7) 



2 



where r is short for 

This substitution does not change the partial sums of series 



«=0 



(8) 



(9) 



for any order A^. If we, however, re-expand these partial sums in powers of g at fixed r up to order A^, and substitute at 
the end r by {co^ - 0.^)/g, we obtain new partial sums 



(10) 



In contrast to E'^, these do depend on the variational parameter Q.. For higher and higher orders, the ^-dependence 
has an increasing valley where the dependence is very weak. It can be found analytically by setting the first derivative 
equal to zero, or, if this equation has no solution, by setting the second derivative equal to zero. One may view this as 
a manifestation of a principle of minimal sensitivity lilOil . The plots are shown in Fig. |2]for odd and even A^. 
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Figure 2: Typical n-dependence of Mh approximations Wjv at 7 = for increasing orders N. The coupling constant has the value g/4 = 0.1. The 
dashed horizontal line indicates the exact energy. 

Even to lowest order, the result is surprisingly accurate. For N = I, the energy we has the linear dependence 

'1 



After the replacement Q and the reexpansion up to power g at fixed r we find 



1 



(11) 



(12) 



14 ■ 4Q ■ 16Q4^ 

In the strong-coupling limit, the minimum lies at Q » c{g/4y^^ where c is some constant and the energy behaves like 



(13) 



The minimum Hes at c = 6^'^ where W' ~ (^/4)'^^ (3/4)"^^^ ^ (g/4)'^^ x 0.681420. The treatment can easily be 
extended to 40 digits 1 11] starting out like (g/4)'/^ x 0.667 986 259 ... . 

The result is shown in for g/4 = 0.1 in Fig. [3] If we plot the minimum as a function of g we obtain the curve 
shown in Fig. [3] The curve has the asymptotic behavior (g/4)'^^ x 0.68142. This grows with the exact power of g 
and has a coefficient that differs only slightly from the accurate value 0.667 986 259 . . . found by other approximation 
procedures il2ll . 



The convergence of the approximations is exponential as was shown in Refs. 11131 Il4 [15f| using the technique 
of order-dependent mapping il7ll . If the asymptotic behavior of E'^ig) and its variational approximation W'^ig) are 
parametrized by 

W'^ig) = ^^{^70 + ^1^^^ +^2^"^ +...}, (14) 
the coefficients bo and bi converge with A^ as shown in Fig. H) The approach is oscillatory (see Fig. |5ll. 
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Figure 4: Asymptotic coefficients and bi of W'^ as a function of the order N. 



3. Quantum Field Theory and Critical Behavior 

When trying to apply the same procedure to quantum field theory, the above procedure needs some important 
modification caused by the fact that the scaling dimensions of fields are no longer equal to the naive dimensions but 
anomalous. This causes the principle of minimal sensitivity to fail |16]. The adaption of the variational procedure 



was done in the textbook lllSll . Let us briefly summarize it using an important class of field theories. 



The energy is an 0(«)-symmetric coupling functional of a n-component field (po in D dimensions 

E[4,Q ] = J d'^xl^ [dM^)f + Y -^oW + |y [M^?f^ , (15) 

where the parameters depend on the distance of the temperature from the critical value Tc'. 

ml = 0{(T-T,y), go^O{(T-T,f) 
The important critical behavior is seen in the correlation function which have the limiting form 



g -|x~x'|/f(r) 



(<^,(x)^,(x'))~— — (16) 
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Figure 5: Oscillations of the strong-coupling coefficient bo. 



where rj is the anomalous field dimension, and ^ is the coherence length which diverges near Tc like ^(T) ~ {T - Tc) 



3.1. Critical Behavior in D — e Dimensions 

The field fluctuations cause divergencies which can be removed by a renormaUzation of field, mass and coupling 
constant to 0, m, and g. This is most elegantly done by assuming the dimension of spacetime to be D = 4 - e, in which 
case the renormalization factor are 



go 
ml 

4 



Z,„{g,e)Z^{g,e) 'ot\ 
Z4,{g, e) ip^. 



The factors have weak-coupling expansions: 



1 + 



n + , 



3e ' 
n + 2 
36e ' 



(n + 8)2 5n + 22 



9e2 



9e 



(17) 
(18) 
(19) 



(20) 
(21) 



n + 2 

Zmig, e) = 1 + — — g + 
3e 



(n + 2)(n + 5) n + 2\ 2 



9e2 



6e 



The dependence of these on the scale parameter ju defines the renormalization group functions 



fi(g,e) = // ^ 
d/u 



7m(g) 



fi dm 
m dfi 



= -e|^ln[gZj(g,e)Z^(g,e)-2]| , 
= -^|-ln[z,„(g,e)Z,(^,e)-'], 



yig) — — — 

(p dfi 



Pig^e) d 

^^lnZ,(g,e). 



(22) 
(23) 
(24) 



At the phase transition go goes to the strong-coupling limit go — > 00. In this limit the renormalized coupling g tends 
to a constant g*, called the fixed point of the theory. 

From the renormalization group functions in the strong-coupling Umit one finds the physical observables at the 
critical point 



^ . *N n + 2 J 
rj = 2y{g*) = — + 



2(n + 8)2 



1 



2[l -y.n{g*)] 2 4(n + 8) 
, , , 3(3n-Hl4) 2 

(n + 8)^ 



1 n + 2 (« + 2)(« + 3)(« -H 20) , 
+ -r, 7:7 e e + 



8(n + 8)3 



(25) 
(26) 
(27) 



The quantity e is the so-called anomalous dimension of the field (pix). 

The e-expansions are divergent and are typically evaluated at the physical value e - I where D = 3 by various 
resummation procedures lll9ll . 

In variational perturbation theory the procedure is different. One rewrites the power series of Eq. ( fTTT i as of the 
renormalized coupling go: 



go ■ 



n + 1 



■^0 + 



(n + 8)2 9n + 42 



9e2 



18e 



gl^ 



(28) 



For the dependence of the renormalized mass on the bare coupling one finds from Eq. ( fTST i 



m^(go) 
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Figure 6: Strong-coupling values of the renormalization group functions for n = I (the so-called Ising universality class). 



and for the anomalous dimension from Eq. ( fT9] l. (l24l l. and dZSl l: 

n + 2 2 (« + 2)(n + 8) 8\ , 
riiSo) = ^1 - -jgo + ■ • ■ • (30) 

Due to the anomalous dimension 77 0, the dependence of the approximations on the variational parameter 
develops no longer a horizontal flat valley (see Appendix A). Instead, the valley turns out to have a slope which can 
only be removed by introducing another parameter q in to substitution rule (|7]). We rewrite the series in ^ as a series 
in g/K^, and replace k by 

yjK^ + (/c2 - K'^) = yjK^ + gr, (31) 

by 

r^(K^- K^)/g. (32) 

As before we re-expand the partial sums of the series in powers of g at fixed r up to power g^ to obtain . After 
this we set /<• — > 1 and plot W'^ as a function of K. By varying q we can make the valley of minimal A'-dependence 
horizontal lfl6ll . 

The asymptotic behavior of the variational parameter K(go) and the critical exponent as a function of go, called 
generically f{go), is now in general 

K(go) = g'"'[co+c,g,'"' + C2g,'"' + ...} 

f(go) = g'"'^{bo + b,g,'"' + b2g,"' + ...], (33) 



In the proof of the exponentially fast convergence in Refs. lilQIil. it was shown that the approach of the 
correct result proceeds as a function of the highest order L of the partial sum as e ''. 
In this way we find from ( l28l l the strong-coupling behavior 1I20I1 



?o) = g* +big^/ +. 



The exponent to is the famous Wegner exponent 112111 . Further we find from 

m\go) , -ir,; ^ 

—J- = ^"^0 +■•■' 



(34) 



(35) 



6 



where the parameter a> and y*, are found from the strong-coupUng Hmits 



e 



-1 -^0 



?"(^o) 



d\n m?- {go) I ml 



d\ngQ 



(36) 



This parameter determines also the divergence of the coherence length in the critical behavior ^(T) ~ (T - Tc) ^: 

y= 1/(2-7;). (37) 

The results are 

1 + ^— e 

2(„.8) ^33^ 



3(3«+14)i; 



2a/i+ („.8): 

They are plotted in Fig. |6]as a function of e 
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Instead of an expansion in D = 4 - e dimensions on may also treat expansions obtained by Nickel 112211 directly in 
D = 3 dimensions. 

3.2. Three-Diinensional Treatment 

If one plots the strong-coupling limits of the series obtained from the partial sums of order L as a function of 
x{L) - e"'^^' " to account for the theoretical approach to the asymptotic limit, one finds for various n ll23i1 : 
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Figure 7: Strong-coupling values for tlie critical exponent v ' (x) as a function of x(L) = e '^^ 

For the critical exponent a characterizing the behavior of the specific heat C ~ |r - Td " of superfluid helium 



near the critical temperature Tc, the strong-coupling limit is OlSll . 



a ^2-3x0.6712 a; -0.0136. (39) 

If we extrapolate the asymptotic behavior expansion coefficients of v up to the 9th order according using the theoret- 
ically known large-order behavior this result can be improved \a a ■x. -0.0129 (see Fig. |8]l. This value agrees 
perfectly with the space shuttle value [gj a - -0.01285 + 0.00038. The experimental result extracted from Fig. [T]and 
the various theoretical numbers obtained from the divergent perturbation series for a are summarized in Fig. |9] 
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Figure 8: Strong-coupling limits of a as a function of jc = e " for 
well with the satellite experiment [Sj]. 
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and 9th order in perturbation theory. The latter limit a » -0.0129 agrees 
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Figure 9: Survey of expeiimental and theoretical values for a. The latter come from resummed perturbation expansions of 0*-theory in 4 - e 
dimensions, in three dimensions, and from high-temperature expansions of XY-models on a lattice. The sources are indicated below. 



4. Shift of the Critical Temperature in Bose-Einstein Condensate by Repulsive Interaction 

A free Bose gas condenses at a critical temperature 



7^(0) ^ 

M 



In 



(40) 



where n is the particle density. A small relative shift of Tc with respect to jf' can be calculated from the general 
formula 

^T. 2 An ^^^^ 



n(0) 



3 n(0) ' 



where rP^ is the particle density in the free condensate and An its change at Tc caused by a small repulsive point 
interaction parametrized by an i-wave scattering length a. For small a, this behaves like lE5lE6l] 



AT, 



= cian^'^ + [C2ln(fln"^) + C2]aV' + O(a'n). 



1/3^ 



,2„2/3 



n(0) 



(42) 



where = -647r^(l/2)/3^(3/2)^^^ ^ 19.7518 can be calculated perturbatively, whereas ci and C2 require non- 
perturbative techniques since infrared divergences at T, make them basically strong-coupling results. The standard 
technique to reach this regime is based on a resummation of perturbation expansions using the renormalization group 
1 13, Ell, first applied in this context by Ref. Q. 
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Using quantum field theory, the temperature shift can be found from the formula 



2 Mr, 



(0) 



AT, 



corresponding in Eq. (l42l l to 



3 n '\ u J" 3 [^(3/2)]4/3 



,1/3 



ci * -1103 



(43) 



(44) 



A calculation of the Feynman diagrams in Fig. [10] yields the following five-loop perturbation expansion for the 
expectation value {(p^/u) |22, 2^ 

Nm N(2 + N)u N (16 + 10 N + N^) , u \^ 

F(u) = «2^T— T— 5 + «3 ,^ 1 — 1 



108 (47r)' 



fl4I 



N{2 + Nf N {4Q + 32N + ^N^+N^) N [44 + 32N + 5 N^) 



324 (4;r)' 



+ 042 



648 {47Ty 



+ 044- 



N(2 + Nf N {U + 32N + 5N^) 



324 (471) 



— + 045- 



324 ot3 (4ny 



+ (243 - 



(-)' 



324 (4ny 



where 02 = log(4/3)/2 ^ 0.143841 and the other constants are only known numerically Bill : 

03 = 0.642144, fl4i = -0.115069, 042 = 3.128107, 043 = 1.63, 044 = -0.624638, 045 = 2.39. 



(45) 



(46) 



Writing the above expansion up to the Lth term as Fl(u) - Ylj^^ yfi(u 1 4nm)' , the expansion coefficients for the relevant 
number of components = 2 are 113 ill : 



/-I = -126.651 X fa = 0, /i = -4.04837 x /z = 2.39701 x 10"\ fi 



-1.80 X 10" 



(47) 



We need the value of the series Fi{u) in the critical limit m — » 0, which is obviously equivalent to the strong-coupUng 
limit of Fiiu). As mentioned above, this limit should be most accurately found with the help of variational perturbation 
theory Slillil. 

If the series were of quantum mechanical origin, we could have found this limit by applying the square-root trick 
O of Ref. |@]. In the present situation where we are only interested in the extreme strong-coupling limit, we would 
form the sequence of truncated expansions Fl(u) for 1, 2, 3 and replace each term 



;/2 



(ulm)' K\\ - \]2 

where the symbol [1 - 1][ is defined as the binomial expansion of (1 - 1)'^ truncated after the Idh term 



,->,; = t(;)(-')' = (-i)'(^;'). 



(48) 



(49) 




Table 1 : Trial functions for the naive quantum-mechanical variational perturbation expansion 
= -0.0596831A'-' - 0.0000322159 A", 
W?^ = -0.0497359A'-' - 0.0000483239 K + 1.51792 lO"*" K^, 
Wf'^ = -0.0435189a:"' - 0.0000604049 K + 3.03584 10"* - .908 10"' . 



The resulting expressions must be optimized in the variational parameter K. They are listed in Table [U The approxi- 
mants WP^3 have extrema ^^^^3"' ^ -0.00277, +0.00405, -0.0029, corresponding, via (gUi, to ci ^ 3.059, -4.46, 3.01. 
These values have previously been obtained in Ref. [29] in a much more complicated way via a so-called cJ-expansion. 
Note the negative sign of the second approximation arising from the fact that an extremum exists only at negative K. 
According to our rules of variational perturbation theory one should, in this case, use the saddle point at positive K 
which would yield - -0.00153 corresponding to c\ » 1.69 rather than -4.46, leading to the more reasonable 
approximation sequence c\ x 3.059, 1.69, 3.01, which shows no sign of convergence. In W^^, there is also a pair of 



complex extrema from which the authors of Ref. II29I1 extract the real part Re W'^^'^jj,piex ~ -0.00134 corresponding to 
ci » 1.48, which they state as their final result. There is, however, no acceptable theoretical justification for such a 
choice 

This lack of convergence is not astonishing since we are dealing with field theory, where the dimensions are 
anomalous and the naive principle of minimal sensitivity breaks down (contrary to ubiquitous statements in the liter- 
ature |34]). The valley in the dependence on the variational parameter is no longer horizontal 1 16]. 

The correct procedure goes as follows: We form the logarithmic derivative of the expansion ( |45] ): 

„„) = « = . 2f (^f . 3f (4f - . . . . ,50) 



1 ^ 



In order for F(u) to go to a constant in the critical limit m — > 0, this function must go to zero in the strong-coupling 
limit M — > 00. Writing the expansion as Piiu) - -I + bi(u/4-nm)', the coeflicients are 

^72 = 0.0639293, ^73 = -0.056778, ^74 = 0.0548799. (51) 

The sums f^iiu) have to be evaluated for m — > 00 allowing for the universal anomalous dimension w by which the 



physical observables of -theories approach the scaling limit Il27lll8r . The approach to the critical point A + B{m/u)'^ 
where a)' = oj/(l — t]/2) Issll . The exponent 77 is the small anomalous dimension of the field while co again the Wegner 
exponent f^T] of renormalization group theory A = a>v. Here it appears in the variational expression for the strong- 
coupling limit which is found Iii32i,.33il by replacing (u/m)' by K'[1 - 1]^^'/^, where q = 2/a>'. Thus we obtain the 
variational expressions 

J-i / J-i 



The first has a vanishing extremum at co'^ = 0.592, the second has neither an extremum nor a saddle point. However, 
a complex pair of extrema lies reasonably close to the real axis at = 0.635 + 0.1 16, whose real part is not far from 
the true exponent of approach a>'^ a; 0.81 ||27, 18|, to which lo'^^ will converge for order L ^ 00 [32]. Given these 
a»'-values, we now form the variational expressions Wl from Fi by the replacement {ulmf — > ^''[l - 1]^^'/^, which 
are 

W2 = /_i(i-^^ + i^2j^-'+/i/:, (54) 



^3 = /-ill-^^ + 



l3-4-'^'-i^1^"^^'(l^f)^^-^^^'' ^''^ 

^4 = /-i(i-g? + g?2-^^^ + ^?4j/:-'+/iJi + ^^ + i^2j^ + /2(i+^)/:2 + /3/:l (56) 
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The lowest function W2 is optimized with the naive growth parameter q - I since to this order no anomalous 
value can be determined from the zero of the beta function (ISOl l. The optimal result is Wj*" - - -\J\og[4-/3]/6/8n^ ~ 
-0.00277 corresponding to ci = 3.06. The next function W3 is optimized with the above determined qj, - and 
yields Wg"" ~ -0.000976 corresponding to ci = 1 .078. Although co'^ is not real we shall insert its real part into W4 and 
find W^"" = -0.000957 corresponding to ci = 1.057. The three values of ci for L = L - 1 - 1,2,3 can well be fitted 

by a function ci a; 1.053 + (see Fig. fTTT l. Such a fit is suggested by the general large-L behavior a + be^'^^' " 
which was derived in Refs. |6]. Due to the smallness of 1 - w' 0.2, this can be replaced by ^ a' + b' /L\ 

Alternatively, we may optimize the functions 1^1,23 using the known precise value of qoo - 2/ll)'^ ^ 2/0.81. Then 
W2 turns out to have no optimum, whereas the others yield W"^' ^ -0.000554, -0.000735, corresponding via Eq. (l44l l 
to ci = 0.580, 0.773. If these two values are fitted by the same inverse power of L, we find ci x 0.83 - 14/U'. From 
the extrapolations to infinite order we estimate ci_oo ~ 0.92 + 0.13. 
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Figure 11: The three approximants for ci plotted against the order of variational approximation L = L - I 
infinite-order limit. 



1,2,3, and extrapolation to the 



This result is to be compared with latest Monte Carlo data which estimate ci k 1.32 ± 0.02 l|36, 37 1. Previous 
theoretical estimates are ci » 2.90 13811 . 2.33 from a 1 /A^-expansion i3% ). 1.71 from a next-to-leading order in a 
expansion [40], 3.059from an inapplicable 5-expansion |41] to three loops, and 1.48 from the same 5-expansion to five 
loops, with a questionable evaluation at a complex extremum [|29.1 and some wrong expansion coeflicients (see 113 ill ). 
Remarkably, our result lies close to the average between the latest and the first Monte Carlo result ci x 0.34 + 0.03 in 
Ref. (4^]. 

As a cross check of the reliability of our theory consider the result in the limit — > 00. Here we must drop the 
first term in the expansion ( |45]) which vanishes at the critical point (but would diverge for — > 00 at finite m). The 
remaining expansion coefficients of /N in powers of Nu/4mn are 



/i 



-6.3591710" 



/2 = 4.7315 10"^ /3 



-3.84146 10" 



(57) 



Using the — > 00 limit of a>' which is equal to 1 implying q - 2m Eqs. ( |55] ) and ( l56b . we obtain the two variational 
approximations 



= -0.00127183/: -H 0.00047315/:^ = -0.00190775/: + 0.00141945/:^ - 0.000384146/:^ 



(58) 



whose optima yield the approximations ci a; 1.886 and 2.017, converging rapidly towards the exact large-A^ result 
2.33 of Ref. IM], with a 10% error. 

Numerically, the first two 1 /A^-corrections found from a fit to large-A^ results obtained by using the known large-A^ 
expression for w' = 1 - 8{8/37t^N) + 2(104/3 - 9;rV2)(8 /37T^Nf \4M produce a finite-A^ correction factor (1 - 3. 1 /A^+ 
30.3/A^2 + ... )^ to be compared with (1 - 0.521 IN + . . .) obtained in Ref. iM . 

Since the large-A^ results can only be obtained so well without the use of the first term we repeat the evaluations 
of the series at the physical value N - 2 without the first term, where the variational expressions for / are 



W2 = /l(l + |)/: + /2/:^ 



W3 = fi\l + '^-q+^-q^\K + f2{l+q)K^+hK^ 



(59) 
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The lowest order optimum lies now at Wj''' — -ff(2 + ^)^/16/|, yielding ci s 0.942 for the exact q - 2/0.81. To 
next order, an optimal turning point of W3 yields ci ^ 1.038. 

At this order, we can derive a variational expression for the determination of oj' using the analog of Eq. ( l50b which 
reads 



_ dlog F(u) f2 u 

P{u) = — = 1 + — - + 

ologM /i m 



■/i /fj 



After the replacement (l48l l we find 



„ f2(l+q/2) 
' /i 

whose vanishing extremum determines co' - 2/q as being 



+ . . . 



.^, = (2^2/i/3//2-l-l) * 0.675, 



(60) 



(61) 



(62) 



leading to ci a: 1.238 from an optimal turning point of Wj. There are now too few points to perform an extrapolation 
to infinite order From the average of the two highest-order results we obtain our final estimate: ci » 1.14±0.11, such 
that the critical temperature shift is 



^^(1.14±0.11)fl„ 



1/3 



(63) 



This lies reasonably close to the Monte Carlo number ci ^ 1 .32 + 0.02. 



5. Membrane Between Walls 

As another example consider a tension-free membrane of bending stiffness between hard walls |44] (see Fig.fTSli. 




Figure 12: Membrane fluctuating between walls with distance d. 



Its thermal fluctuations are described by a functional integral over a Boltzmann factor 

V/ J-d/2 

where h(x) is the height function of the membrane and E is the bending energy 

n2 



(64) 



(65) 



This functional integral has not been solved exactly, in spite of its simplicity. It can, however, be approximated by the 
functional integral 



c£)hg-[E+V(x)]lhT 
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(66) 




m ^ ■ 



-dl2 h d/2 

Figure 13: Softened hard-wall potential which becomes infinitely hard in the Umit m —> 



in which the height fluctuates between -co and oo in a potential (see Fig. [TsT i 



y(jc) = m4— tan^ — 



(67) 



This problem can be solved perturbatively yielding Z - e , where A is the area of the membrane and / has, to order 
the series 



r 



1 + + 



8 m^dP- 64 



The hard-wall limit m amounts to the strong-coupling limit of this series. 
We expand the potential © into a power series 

V(h) = m^— +m^ — \-h* + + + /i'^ + /;"' + 

2 c/2 1 3 90 315 c/4 14175 ^6 467775 



If we denote the interaction terms by 



4 ^i. 



(68) 



(69) 



(70) 



and calculate the Feynman diagrams shown in Fig. [14] The functional integral ( l64b can be expressed as an exponential 



1 



1 



F'* - - O + 300 + 15 + - 



^2 000 +24^ 



+ 



6 



+105 ^ +- (540 ^ + 360 ) 
1728 15 + 3456 R + 1728 ^ + 2592 OCOO 



Figure 14: Feynman diagrams in the perturbative expansion of the free energy of the Membrane between walls up to the order N = A. 

Z — e where A is the area of the membrane and 



■' 2 



Using the Bender-Wu recursion relations 114611 . we can express the coefficients in terms of sk as 

2 2 4 6 

^'-T^ Z^^4 - (21e^ - ISfie) + (333e^-360e4fi6 + 105e8) 



8 

-Y^^(308854-44880e^e6+6990e^-Hl512e4e8 + 3780ei())- 
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(71) 



The hard-wall result is obtained in the limit m — > 0, which is the strong-coupling limit of the series ( fTTT i. 



6. Variational Perturbation Theory of T\inneling 

None of the presently known resummation schemes ifiol. 18 1 is able to deal with non-Borel-summable series. Such 
series arise in the theoretical description of many important physical phenomena, in particular tunneling processes. 
In the path integral, these are dominated by non-perturbative contributions coming from nontrivial classical solutions 
called critical bubbles ll45[ |6ll or bounces f47ll . and fluctuations around these. 

A non-Borel-summable series can become Borel-summable if the expansion parameter, usually some coupling 
constant g, is continued to negative values. In this way, non-Borel-summable series can be evaluated with any desked 
accuracy by an analytic continuation of variational perturbation theory in the complex g-plane. This impUes 

that variational perturbation theory can give us information on non-perturbative properties of the theory. 

6.1. Test of Variational Perturbation Theory for Simple Model of Non-Borel-summable Expansions 
The partition function Z(g) of the anharmonic oscillator in zero space-time dimensions is 

Z(g)=^ r exp(-xV2 - g //4) dx ^ ^''P^.!^!^V i/4(l/8g) , (72) 

where Ky(z) is the modified Bessel function. For small g, the function Z(g) has a divergent Taylor series expansion, to 
be called weak-coupling expansion: 



Z 

weak 



•W . ^ V ' >K . .,,r(2[+l/2) 

" (g)^ ) ai g', with ai = (-1)' — — — . (73) 



For g < 0, this is non-Borel-summable. For large |g| there exists a convergent strong-coupling expansion: 



L 



Z^.,(8) = 8-"' Z b, g-"\ with b, = (-1)' (74) 



/=o 



2/! 



As is obvious from the integral representation (l72l l. Z(g) obeys the second-order diflerential equation 

l6g^Z"(g) + 4(1 + 8g)Z'(g) + 3Z(g) = 0, (75) 

which has two independent solutions. One of them is Z(g), which is finite for g > with Z(0) = aq- The weak- 
coupling coefficients ai in ( |73] | can be obtained by inserting into (iTST i the Taylor series and comparing coefficients. 
The result is the recursion relation 

16/(/+l) + 3 

fl;+l — ai. (7o) 

A similar recursion relation can be derived for the strong-coupling coefficients bi in Eq. ( l74l l. We observe that 
the two independent solutions Z(g) of dTST l behave like Z(g) oc g" for g — > oo with the powers a = -1/4 and -3/4. 
The function ( 172] ) has a = - 1 /4. It is convenient to remove the leading power from Z(g) and define a function ^(x) 
such that Z(g) = g"'^"* ^(g"'^^). The Taylor coefficients of ^(x) are the strong-coupling coefficients bi in Eq. ( l74b . The 
function ^(x) satisfies the differential equation and initial conditions: 

4C'(x) - 2x('(x) - ((x) = 0, with ^(0) = bo and ^'(0) = bi. (77) 

The Taylor coefficients bi of ^(x) satisfy the recursion relation 
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Analytic continuation of Z(g) around g = oo to the left-hand cut gives: 



1=0 



-(21+1) 



for g > 0, 



(79) 
(80) 



so that we find an imaginary part 



ImZ(-g) = -(4g)-'/4^Z7,(-^)-'/2sin 

1=0 



-'-^(21+ I) 



1=0 



where 



It is easy to show that 



so that 



21+1 



1+2 



m + + 2) 



Pi 



^A^' = f(x)exp(-^V4), 



1 °° 
ImZi-g) = -^g'"' exp(-l/4g) ^ big 



V2' ,=„ 
From this we may re-obtain the weak-coupling coefficients a; by means of the dispersion relation 



Z(g) 



1 r°°ImZ(-z) 



^ Jo z + , 



dz 



exp(-l/4z)z->/2-i/4 



1 y ^ r exp(-l/4z 
ny/2 Jo z + 



dz. 



(81) 
(82) 

(83) 
(84) 

(85) 

(86) 
(87) 



Indeed, replacing l/(z + g)hy exp (-x(z + g)) dx, and expanding exp (-x g) into a power series, all integrals can 
be evaluated to yield: 



^ CO CO 

z(g) ^ 2'bj ^ (-gym + j/2 + 1/4) 



(88) 



j=0 1=0 

Thus we find for the weak-coupling coefficients a/ an expansion in terms of the strong-coupling coefficients 

(-1)' 

Inserting bj from Eq. ( |74] |. this becomes 



ai 



y 2>bjY(l + jl2 + \IA). 



(89) 



ai - 



coinciding with ( l73T l. 



t%f^^JSdl nji2 . ./4)ra . ;,2 . 1/4) = (-.). IHli/Hl , 
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(90) 



Let us now apply variational perturbation theory to the weak-coupling expansion ( l73T l. We have seen in Eq. ( |79] l, 
that the strong-coupling expansion can easily be continued analytically to negative g. This continuation can, however, 
be used for an evaluation only for sufficiently large \g\ where the strong-coupling expansion converges. In the tunneling 
regime near the tip of the left-hand cut, the expansion diverges. Let us show that an evaluation of the weak-coupling 
expansion according to the rules of variational perturbation theory continued into the complex plane gives extremely 
good results on the entke left-hand cut with a fast convergence even near the tip at g - Q. 

The Lth variational approximation to Z(g) is given by (see EllllllllTl]) 

zfj(g,Q) = Q'' e/cr), (91) 



with 

where q - 2/a) - A, p - -\ and 



cr = 0«-2(Q2_i)/^, (92) 



1=0 ^ ' 

In order to find a valley of minimal sensitivity, the zeros of the derivative of z['^l(g, Q) with respect to Q are needed. 
They are given by the zeros of the polynomials in cr; 

P'^'icr) = g a,(p -lq + 2l- 2L)(^^ -^^J/^)(-cr)^-' = 0, (94) 



since it can be shown 013111511 that the derivative depends only on cr: 



dZ^l(g,£l) 



n"-' (^^)' P^^\(r) . (95) 



Consider in more detail the lowest non-trivial order with L - 1. From Eq. ( |94l l we obtain 




Figure 15: Plot of the 1st- and 2nd-order calculation for the non-Borel-summable region of g < 0, where the function has a cut with non-vanishing 
imaginary part: imaginary (left) and real parts (right) of zi,!^.(g) (dashed curve) and Z*^(g) (solid curve) are plotted against g and compared with 
the exact values of the partition function (dotted curve). The root of i92i giving the optimal variational parameter Q. has been chosen to reproduce 
the weak-coupling result near g = 0. 



<T corresponding to Q. - - + ^^s) ■ (96) 

In order to ensure that our method reproduces the weak-coupling result for small g, we have to take the positive sign in 
front of the square root. In Fig.[T5]we have plotted Zvar(g) (dashed curve) and Z^^^ig) (solid curve) and compared these 
with the exact result (doted curve) in the tunneling regime. The agreement is quite good even at these low orders llsill . 
Next we study the behavior of zl^'^lig) to higher orders L. For selected coupling values in the non-Borel-summable 
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region, g — -.01, -.1, -1, -10, we want to see the error as a function of the order. We want to find from this model 
system the rule for selecting systematically the best zero of P*^'(cr) solving Eq. ( |94] l, which leads to the optimal value 
of the variational parameter Q. For this purpose we plot the variational results of all zeros. This is shown in Fig.fTSl 
where the logarithm of the deviations from the exact value is plotted against the order L. The outcome of different 
zeros cluster strongly near the best value. Therefore, choosing any zero out of the middle of the cluster is reasonable, 
in particular, because it does not depend on the knowledge of the exact solution, so that this rule may be taken over to 
realistic cases. 
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Figure 16: Logarithm of deviation of the variational results from exact values log - ZexactI plotted against the order L for different g < in the 
non-Borel-summable region. All complex optimal il's have been used. 



A(L) 



Figure 17: Logarithm of deviation of variational results from exactly known value A(L) = log \Zy.J - Z^xaaV plotted against the order L for g = 10 
in Borel-summable region. The real positive optimal n have been used. There is only one real zero of the first derivative in every odd order L and 
none for even orders. There is excellent convergence A(L) ^ 0.02 exp (-0.73L) for L — » cxj. 



We wish to emphasize, that for the Borel-summable domain with g > 0, variational perturbation theory has the 
usual fast convergence in this model. In fact, for g - 10, probing deeply into the strong-coupling domain, we find 
rapid convergence like A(L) ^ 0.02 exp (-0.73L) for L — > oo, where A(L) — log |Zvai - ZexactI is the logarithmic error 
as a function of the order L. This is shown in Fig. [T7] Furthermore, the strong-coupling coefficients hi of Eq. ( l74l l 
are reproduced quite satisfactorily. Having solved P^^\cr) - for cr, we obtain Q*^'(g) by solving Eq. ( |92l ). Inserting 
this and ( |93] ) into ( 1911 1. we bring g'^* Zvar(g) into a form suitable for expansion in powers of The expansion 

coefficients are the strong-coupling coefficients b^j'^ to order L. In Fig. [18] we have plotted the logarithms of their 
absolute and relative errors over the order L, and find very good convergence, showing that variational perturbation 
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Figure 18: Relative logarithmic error A, = log |1 - /b^ \ on the left, and the absolute logarithmic error A„ = log 
plotted for some strong-coupling coefficients bi with I = 0,4, 8, 12, 16, 20 against the order L. 



, (exact) 



I on the right. 



theory works well for our test-model Z(g). 

A better selection of the optimal Q values comes from the following observation. The imaginary parts of the 
approximations near the singularity at ^ = show tiny oscillations. The exact imaginary part is known to decrease 
extremely fast, like exp (l/4g), for g — > 0-, practically without oscillations. We can make the tiny oscillations more 
visible by taking this exponential factor out of the imaginary part. This is done in Fig. [19] The oscillations differ 
strongly for different choices of Q*^* from the central region of the cluster To each order L we see that one of them 
is smoothest in the sense that the approximation approaches the singularity most closely before oscillations begin. If 
this Q*^* is chosen as the optimal one, we obtain excellent results for the entire non-Borel-summable region g < 0. 
As an example, we pick the best zero for the L = 16th order Fig.[T9]shows the normalized imaginary part calculated 




Figure 19: Normalized imaginary part Im[z'!^'(g) exp (-l/4g)] as a function of g based on six different complex zeros (thin curves). The fat curve 
represents the exact value, which is Zexactfe) - -0.7071 + .524g - 1.78g^. Oscillations of varying strength can be observed near g = 0. Curves 
A and C carry most smoothly near up to the origin. Evaluation based on either of them yields equally good results. We have selected the zero 
belonging to curve C as our best choice to this order L = 16. 



to this order, but based on different zeros from the central cluster Curve C appears optimal. Therefore we select the 
underlying zero as our best choice at order L = 16 and calculate with it real and imaginary part for the non-Borel- 
summable region -2 < g < -.008, to be compared with the exact values. Both are shown in Fig.|20l where we have 
again renormalized the imaginary part by the exponential factor exp(-l/4g). The agreement with the exact result 
(solid curve) is excellent as was to be expected because of the fast convergence observed in Fig. [16] It is indeed 
much better than the strong-coupling expansion to the same order, shown as a dashed curve. This is the essential 



improvement of our present theory as compared to previously known methods probing into the tunneling regime 115 111 . 
This non-Borel-summable regime will now be investigated for the quantum-mechanical anharmonic oscillator 
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-2 -4 log(-s) -2 -4 log(-s) 

Figure 20: Normalized imaginary part Im[Zi,|,j'(g) exp (-l/4g)] to the left and the real part Re[Zy!,f to the right, based on the best zero C from 
Fig. [19] are plotted against log \g\ as dots. The solid curve represents the exact function. The dashed curve is the 16th order of the strong-coupling 
expansion Z^j^* (g) of equation )74t . 



6.2. Tunneling Regime of Quantum-Mechanical Anharmonic Oscillator 

The divergent weak-coupling perturbation expansion for the ground state energy of the anharmonic oscillator in 
the potential V{x) - /2 + g to order L 

L 

Ceak(^) = Z«'^'' ^97) 

/=0 

wherea/ = (1/2, 3/4, -21/8, 333/16, -30885/128, . . . ), is non-Borel-summable for ^ < 0. It may be treated in the 
same way as Z(g) of the previous model, making use as before of Eqs. (I9ni-(l94]i. provided we set p = 1 and o) -Ijli, 
so that q - li, accounting for the correct power behavior Eo{g) cc g^/^ for g — > oo. According to the principle of 
minimal dependence and oscillations, we pick a best zero for the order L = 64 from the cluster of zeros of Piicr), and 
use it to calculate the logarithm of the normaUzed imaginary part: 

fig) : = log f Ef^^ig)] -l/3g. (98) 

This quantity is plotted in Fig.[2T]against \og(-g) close to the tip of the left-hand cut for -.2 < g < -.006. Comparing 
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Figure 21: Logarithm of the imaginary part of the ground state energy of the anharmonic oscillator with the essential singularity factored out for 
better visualization, 1(g) = log |^ ^J-ng/l - l/3g, plotted against small negative values of the coupling constant -0.2 < g < -.006 where 

the series is non-Borel-summable. The thin curve represents the divergent expansion around a critical bubble of Ref. | 52]. The fat curve is the 22nd 
order approximation of the strong-coupling expansion, analytically continued to negative g in the sliding regime calculated in Chapter 17 of the 
textbook lij. 

our result to older values from semi-classical calculations lf52|] 

fig) = big - b2g^ + b3g' - b.g" + ..., (99) 

with 

bi = 3.95833 b2 = 19.344 b^ = 174.21 b^ = 2177 , (100) 
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-2 -3 -4 -5 log(-g) 

Figure 22: Logarithm of tlie normalized imaginary part of the ground state energy log ( ^f-TrgJl EQ^^.^(g)) - 1 /3g, plotted against log (-g) for orders 
L = 4, 8, 16, 32 (curves). It is compared with the corresponding results for L = 64 (points). This is shown for small negative values of the coupling 
constant -0.2 < g < -.006, i.e. in the non-Borel-summable critical-bubble region. Fast convergence is easily recognized. Lower orders oscillate 
more heavily. Increasing orders allow closer approach to the singularity aX g = 0-. 



shown in Fig. |2T] as a thin curve, we find very good agreement. This expansion contains the information on the 
fluctuations around the critical bubble. It is divergent and non-Borel-summable for g < 0. In Appendix B we have 
rederived it in a novel way which allowed us to extend and improve it considerably. 

Remarkably, our theory allows us to retrieve the first three terms of this expansion from the perturbation expansion. 
Since our result provides us with a regular approximation to the essential singularity, the fitting procedure depends 
somewhat on the interval over which we fit our curve by a power series. A compromise between a sufficiently long 
interval and the runaway of the divergent critical-bubble expansion is obtained for a lower limit g > -.0229 ± .0003 
and an upper limit g = -0.006. Fitting a polynomial to the data, we extract the following first three coefficients: 

fei = 3.9586 ± .0003 /?2 = 19.4±.12 ^73 = 135 + 18 . (101) 

The agreement of these numbers with those in ( |99] l demonstrates that our method is capable of probing deeply into 
the critical-bubble region of the coupling constant. 

Further evidence for the quality of our theory comes from a comparison with the analytically continued strong- 
coupling result plotted to order L = 22 as a fat curve in Fig. |2T| This expansion was derived by a procedure of 
summing non-Borel-summable series developed in Chapter 17 of the textbook |6]. It was based on a two-step process: 
the derivation of a strong-coupling expansion of the type (l74l i from the divergent weak-coupling expansion, and an 
analytic continuation of the strong-coupling expansion to negative g. This method was applicable only for large 
enough coupling strength where the strong-coupling expansion converges, the so-called sliding regime. It could not 
invade into the tunneling regime at small g governed by critical bubbles, which was treated in |6] by a separate 
variational procedure. The present work fills the missing gap by extending variational perturbation theory to all g 
arbitrarily close to zero, without the need for a separate treatment of the tunneling regime. 

It is interesting to see, how the correct limit is approached as the order L increases. This is shown in Fig. |22] 
based on the optimal zero in each order. For large negative g, even the small orders give excellent results. Close to the 
singularity the scaling factor exp (- 1 /3g) will always win over the perturbation results. It is surprising, however, how 
fantastically close to the singularity we can go. 

6.3. Dynamic Approach to the Critical-Bubble Regime 

Regarding the computational challenges connected with the critical-bubble regime of small g < 0, it is worth to 
develop an independent method to calculate imaginary parts in the tunneling regime. For a quantum-mechanical sys- 
tem with an interaction potential g V{x), such as a the harmonic oscillator, we may study the effect of an infinitesimal 
increase in g upon the system. It induces an infinitesimal unitary transformation of the Hilbert space. The new Hilbert 
space can be made the starting point for the next infinitesimal increase in g. In this way we derive an infinite set of first 
order ordinary differential equations for the change of the energy levels and matrix elements (for details see Appendix 
C): 
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Figure 23: inmy part of the ground state energy of the anharmonic oscillator as solution of the coupled set of differential equations U02K truncated 
at the energy level of n = 64 (points), compared with the corresponding quantity from the L = 64th order of non-Borel-summable variational 
perturbation theory (curve), both shown as functions of the coupling constant g. 



E'„(8) -V„„(8l (102) 



T/' ^ \ V ^'"k{g)Vkn{g) , V ^i"k(8)yknKg) ,,„^. 

E„(g) - Ek(g) E„(g) - Ek(g) 



This system of equations holds for any one-dimensional Schroedinger problem. Individual differences come from the 
initial conditions, which are the energy levels £',,(0) of the unperturbed system and the matrix elements y„m(0) of the 
interaction V(x) in the unperturbed basis. For a numerical integration of the system a truncation is necessary. The 
obvious way is to restrict the Hilbert space to the manifold spanned by the lowest eigenvectors of the unperturbed 
system. For cases like the anharmonic oscillator, which are even, with even perturbation and with only an even state 
to be investigated, we may span the Hilbert space by even basis vectors only. Our initial conditions are thus for 
« = 0, 1, 2, NI2: 

£2„(0) =2n + 1/2 (104) 
V2„,2m =0 if m < or m > N jl (105) 
V2„,2„(0)=3(8«2+4« + l)/4 (106) 
V2„,2„±2(0) =(4n + 3) V(2« + 1)(2« + 2)/2 (107) 

V2,,2„±4(0) = V(2n + l)(2n + 2)(2« + 3)(2« + 4)/4 (108) 

(109) 



For the anharmonic oscillator with a V{x) - potential, all sums in equation ( I102l i are finite with at most four terms 
due to the near-diagonal structure of the perturbation. 

In order to find a solution for some ^ < 0, we first integrate the system from to \g\, then around a semi-circle 
g = \g\ exp (itp) from ip = to ip - n. The imaginary part of Eo{g) obtained from a calculation with = 64 is shown 
in Fig.|23] where it is compared with the variational result for L - 64. The agreement is excellent. It must be noted, 
however, that the necessary truncation of the system of differential equations introduces an error, which cannot be 
made arbitrarily small by increasing the truncation limit A^. The approximations are asymptotic sharing this property 
with the original weak-coupling series. Its divergence is, however, reduced considerably, which is the reason why we 
obtain accurate results for the critical-bubble regime, where the weak-coupling series fails completely to reproduce 
the imaginary part. 

7. Hydrogen Atom in Strong Magnetic Field 

A point particle in D dimensions with a potential V{x) and a vector potential A(x) is described by a Hamiltonian 



i/(p,x) = 



Ire 1^ 
p--A(x) ■ 

c 
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47r|x| 



(110) 



The quantum statistical partition function is given by the euclidean phase space path integral 



with an action 



and the path measure 



J?1[P,X]= fl'T[-/p(r)-x(T)+i/(p(T),x(T))], 

Jo 



(111) 



(112) 



(113) 



The parameter p - l/ksT denotes the usual inverse thermal energy at temperature T, where kg is the Boltzmann 
constant. From Z we obtain the free energy of the system: 



F = --InZ. 

/3 



(114) 



Applying variational perturbation theory to the path integral (II 1 II ) leads to a variational binding energy 1154 1 



defined by s(B) = B /2 - E(B) in atomic natural with /i = 1, M - \,e- 1, energies in units of2Ryd= e'^M^/ti? 



in BO/ 77\ B" InLl 
(B) = 1 + - a/ —h{i]) 



T]\ B^ 



with 



1 



1- yr^ 



:ln 



Here we have introduced variational parameters 

2Q„ 

Q = 

Extremizing the energy with respect to these yields the conditions 



2Q|| 

77 = < 1, 



(115) 



(116) 



(117) 



Q. 1 
2nri 1 - J] 



1 ^j] 5 _^ 1 



1 + 



1 1 



:ln 



1- yr^ 

2 VT^'"" 1 + VT^J 

1 . 1- vr^ 



:ln 



= 0. 



4 8 4Q2 2 \ 2nQ i + 

Expanding the variational parameters into perturbation series of the square magnetic field B^, 

CO CO 



(118) 



(119) 



and inserting these expansions into the self-consistency conditions ( II 181 l and (11181 ) we obtain order by order the 
coefficients given in Table |2] Inserting these values into the expression for the binding energy (II 15l l and expand with 
respect to B^, we obtain the perturbation series 



e^^\B)^^-J]s„B'". 



(120) 



n=0 



The first coeflicients are also given in Table |2l We find thus the important result that the first-order variational 
perturbation solution possesses a perturbative behavior with respect to the square magnetic field strength B^ in the 
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Table 2: Perturbation coefficients up to order 6 for the weak-field expansions of the variational parameters and the binding energy in comparison 
to the exact ones of Ref. (S^l . 



n 







1 


2 


3 


Tin 




1.0 


"^W ~ -0-5576 


16828965T-' , :!rm 
1258815488 ~ ^■J"^-' 


3886999332075J-'' ^ ^^tn 
884272562962432 ~ 


Q„ 


32 


« 1.1318 


^ 1.3885 


I293975T-' n mnoi 


5244316671871= . c on-7-7 


971 


19668992 ~ ^■OJi'O^ 


27633517592576 ~ 




4 

in 


X -0.4244 
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weak- field limit thus yielding the correct asymptotic. The coefficients differ in higher order from the exact ones but 
are improved by variational perturbation theory |6]. 
In a strong magnetic field one has 

» 2Q||, Q|| « B 



and the variational expression simplifies to 



,(1) 
•'"±,"11 



B 



B 



— -H — + ^+ W^ln 



4Qi 



n 20^ j 



which is minimal at 



— -(lnQ|| - InQi + 2 - ln2), 



2a — +Ba 1 +4 



Expanding the second conditions as 



/q7 On 

y n 7tB tt^B^ 



and inserting only the first two terms into the first condition (11231) . we neglect terms of order 1 /B, and find 



Q|| ~ — (inB - lnQ[|" + ln2 - 2) . 
To obtain a tractable approximation for Oy, we perform some iterations starting from 



\n2Be 



-2 



(121) 
(122) 

(123) 
(124) 

(125) 
(126) 

(127) 



Reinserting this on the right-hand side of Eq. (1126b . one obtains the second iteration -jOn . We stop this procedure 
after an additional reinsertion which yields 



/nP' = |ln2Be"2-21n 



^|l„2fi.--21„(^,„26.-)|]). 



(128) 



The reader may convince himself that this iteration procedure indeed converges. For a subsequent systematical ex- 
traction of terms essentially contributing to the binding energy, the expression ( I128l l is not satisfactory. Therefore 
it is better to separate the leading term in the curly brackets and expand the logarithm of the remainder Then this 



procedure is applied to the expression in the square brackets and so on. Neglecting terms of order In ^B, we obtain 

Jflf^ ^ -^(\n2Be-^ + \n- - 21nln 2Be-A. ( 1 29) 

" s/tt \ 4 / 
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Table 3: Example for the competing leading six terms in Eq. 1134) at B = 10^ Bq ~ 2.35 X lO''* G. 

(l/ff)ln^B -(4/ff)ln B Inln B (4/;r)ln"lnB -(4h/n)MnB [2(h + 2)/n]ln B b^/ir 
42.1912 -35.8181 7.6019 4.8173 3.3098 0.7632 



The double-logarithmic term can be expanded in a similar way as described above 

ln2-2 " 



ln]n2Be-^ = In 



InB 1 + 



InB 



ln2-2 l(ln2-2)2 ^ , 

= InlnB + + OilrT^B). 

InB 2 \n^B 



Thus the expression ( |129t may be rewritten as 



2 



2a a 

lnB-21nlnB+ + 

InB 



In^B 



+ b\ + Oiln-^B) 



(130) 



(131) 



with abbreviations 



fl = 2-ln2 a: 1.307, = In- - 2 a -1.548. (132) 

The first observation is that the variational parameter is always much smaller than in the high B-field limit. 



(133) 



Thus we can further simplify the approximation ( 1125b by replacing 



a B 



1 + - A/ — 
B V 7T 



without affecting the following expression for the binding energy. Inserting the solutions ( 1131b and ( 1133b into the 
equation for the binding energy (1122b and expanding the logarithmic term once more as described, we find up to the 
order In^^B: 



(B) = i ( In^B - 4 In B Inln B + 4 In^ln B - 4^ Inln B + 2(^ + 2) In B + - — [8 In^ln B - 8^ Inhi B + 2b^] 
7t\ In B l ij 

+0(ln-^B) (134) 



.(1) 



Note that the prefactor 1 /n of the leading In B-term differs from a value 1 /2 obtained by Landau and Lifschitz [56 
Our different value is a consequence of using a harmonic trial system. The calculation of higher orders in variational 
perturbation theory would improve the value of the prefactor 

At a magnetic field strength B = lO^'Bo, which corresponds to 2.35 x 10'"T = 2.35 x lO'^^G, the contribution from 
the first six terms is 22.87 [2 Ryd]. The next three terms suppressed by a factor In^'B contribute -2.29 [2 Ryd], while 
an estimate for the In^^B-terms yields nearly -0.3 [2 Ryd]. Thus we find 



£^"(10^) = 20.58 ±0.3 [2 Ryd]. 



(135) 



This is in very good agreement with the value 20.60 [2 Ryd] obtained from an accurate numerical treatment fs^. 

Table |3]lists the values of the first six terms of Eq. (1134b . This shows in particular the significance of the second- 
leading term -(4/7r)ln B Inln B, which is of the same order of the leading term (1 /7r)ln^B but with an opposite sign. In 
Fig. |24] we have plotted the expression 

(136) 



edB) = i In^B 



from Landau and Lifschitz [56] to illustrate that it gives far too large binding energies even at very large magnetic 
fields, e.g. at 2000Bo oc lO'^G. 

This strength of magnetic field appears on surfaces of neutron stars (10'" - lO'^G). A recently discovered new 
type of neutron star is the so-called magnetar In these, charged particles such as protons and electrons produced by 
decaying neutrons give rise to the giant magnetic field of 10'^ G. Magnetic fields of white dwarfs reach only up to 
10^ - 10^ G. All these magnetic field strengths are far from realization in experiments. The strongest magnetic fields 
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ever produced in a laboratory were only of the order 10^ G, an order of magnitude larger than the fields in sun spots 
which reach about 0.4 x 10"^ G. Recall, for comparison, that the earth's magnetic field has the small value of 0.6 G. 

The nonleading terms in Eq. (1134b give important contributions to the asymptotic behavior even at such large 
magnetic fields, as we can see in Fig. |24] It is an unusual property of the asymptotic behavior that the absolute 
value of the difference between the Landau-expression (1136b and our approximation ( 1134b diverges with increasing 
magnetic field strengths B, only the relative difference decreases. 




Figure 24: Ground state energy E{B) of hydrogen in a strong magnetic field The dotted figure on the left is Landau's old upper hmit. On the light- 
hand side our curve is compared with the accurate values (dots |57, 58]). It also shows various lower-order approximations within our procedure. 
The quantity e(B) is the binding energy defined by £(B) = B/2 ~ E(B). All quantities are in atomic natural units h = \, M = \, e = \, energies in 
units of 2Ryd= e'*M^/h^. 



8. Appendix A: Modification of Principle of Minimal Sensitivity 

The naive quantum mechanical variational perturbation theory has been used by many authors under the name 
^-expansion. This name stems from the fact that one may write the Hamiltonian of an anharmonic oscillator 

+ + (137) 

2M 2 4 



alternatively as 



2M 2 



(138) 



and expand the eigenvalues systematically in powers of 6. Each partial sum of order L is evaluated at 6 - I and 
extremized in Q. It is obvious that this procedure is equivalent the re-expansion method in Section|2] 

As mentioned in the text and pointed out in [ 16], such an analysis is inapplicable in quantum field theory, where the 
Wegner exponent w is anomalous and must be determined dynamically. Most recently, the false treatment was given to 



the shift of the critical temperature in a Bose-Einstein condensate caused by a small interaction 115011291 14111 . We have 
seen in Section|4]that the perturbation expansion for this quantity is a function of g //i where jj is the chemical potential 
which goes to zero at the critical point, we are faced with a typical strong-coupling problem of critical phenomena. In 
order to justify the application of the 5-expansion to this problem, BR [50] studied the convergence properties of the 
method by applying it to a certain amplitude A{g) of an (9(A^)-symmetric ip*-field theory in the limit of large A^, where 
the model is exactly solvable. 

Their procedure must be criticized in two ways. First, the amplitude A(^) they considered is not a good candidate 



for a resummation by a (5-expansion since it does not possess the characteristic strong-coupling power structure US | 
of quantum mechanics and field theory, which the final resummed expression will always have by construction. The 
power structure is disturbed by additional logarithmic terms. Second, the ^-expansion is, in the example, equivalent 
to choosing, on dimensional grounds, the exponent co - 2 in fisl, which is far from the correct value s; 0.843 to 
be derived below. Thus the 5-expansion is inapplicable, and this explains the problems into which BR run in their 
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resummation attempt. Most importantly, they do not find a well-shaped plateau of the variational expressions A^^\g, z) 
as a function of z which would be necessary for invoking the principle of minimal sensitivity. Instead, they observe 
that the zeros of the first derivatives d~A^^\g,z) run away far into the complex plain. Choosing the complex solutions 
to determine their final resummed value misses the correct one by 3% up to the 35th order. 

One may improve the situation by trying out various dififerent w-values and choosing the best of them yielding an 
acceptable plateau in A(g,z)- This happens for oj x 0.843. However, even for this optimal value, the resummation 
result never converges to the correct limit. For A{g) the error happens to be numerically small, only 0. 1 %, but it will 
be uncontrolled in physical problems where the result is unknown. 

Let us explain these points in more detail. BR consider the weak-coupling series with the reexpansion parameter 

6: 

A{6,g)^-Y{ ■J=\ai, where ai= \ K{x)f'{x)dx, (139) 

Jo 

with 

4x^ 2 X 

7t(1 + x^y X 2 

The geometric series in ( 11391 ) can be summed exactly, and the result may formally be reexpanded into a strong- 
coupling series in h = yl^^/(6 g): 

A(6,g)= I K(x) dx^y b„,{-hT , where b„,^ \ K(x)f-"'(x) dx . (141) 

Jo ^/^^ + 6gfix) tt> 

nOC 

The strong-coupling limit is found for h where A ^ bo - Jq dx K(x) - 1 . The approach to this limit is, however, 
not given by a strong-coupling expansion of the form (1141b . This would only happen if all the integrals b^ were to 
exist which, unfortunately, is not the case since all integrals for b,,, with m > Q diverge at the upper limit, where 

/(x) = - arctan ^ ~ - . (142) 
X 2 X 

The exact behavior of A in the strong-coupling limit /i — » is found by studying the effect of the asymptotic n/x- 
contribution of f(x) to the integral in (I141l i. For f(x) - n/xwe obtain 



r 

Jo 



1 n'^ + 2nh - nh + 2h + Anhloghin 

1 + n/j(x) (n + h) 



The logarithm of h shows a mismatch with the general asymptotic form of the result II15I1 . which and prevents the 
expansion ( 11391 ) to be a candidate for variational perturbation theory. 

We now explain the second criticism. Suppose we ignore the just-demonstrated fundamental obstacle and follow 



the rules of the 5-expansion, defining the Lth order approximant A(6, oo) by expanding ( 1139b in powers of 6 up to order 
5^, setting 6-1, and defining z = g- Then we obtain the Lth variational expression for bo: 



L 



bfic.,z)^Y.a,ji^~l-'!"]^ (144) 

with w = 2, to be optimized in z- This w-value would only be adequate if the approach to the strong-coupling limit 
behaved like A + B/h^ + . . ., rather than ( 1143b . This is the reason why BR find no real regime of minimal sensitivity 
onz. 

Let us attempt to improve the situation by determining o) dynamically by making the plateau in the plots of 
A^^\a>,h) versus h horizontal for several different w-values. The result is w ^ 0.843, quite far from the naive value 
2. This value can also be estimated by inspecting plots of A^^\a), h) versus h for several different w-values in Fig.l25l 
and selecting the one producing minimal sensitivity. It produces reasonable results also in higher orders, as is seen in 
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-0.5 t ^ ^ ^ ^ 

0.5 1 1.5 2.5 

Figure 25: Plot of 1 - b'j^\oj,z) versus z for L = 10 and w = 0.6, 0.843, 1, 2 . The curve with oj = 0.6 shows oscillations. They decrease with 
increasing (o and becomes flat at about <o = 0.843. Further increase of (o tilts the plateau and shows no regime of minimal sensitivity. At the same 
time, the minimum of the curve rises rapidly above the correct value of 1 - i)o = 0, as can be seen from the upper two curves for = 1 and (o = 2, 
respectively. 
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Figure 26: Left-hand column shows plots of 1 - fcj, '(w, z) for L = 10, 17, 24, 31, 38, 45 with &> = 2 of (5-expansion of BR, right-hand column 
with optimal = 0.843. The lower row enlarges the interesting plateau regions of the plots above. Only the right-hand side shows minimal 
sensitivity, and the associated plateau lies closer to the correct value 1 - fco = than the minima in the left column by two orders of magnitude. 
StiU the right-hand curves do not approach the exact limit for L -» oo due to the wrong strong-coupling behavior of the initial function. 
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Figure 27: Deviation of 1 - fe], pjateau^'^ ~ 0.843) ftom zero as a function of the order L. Asymptotically the value -.001 136 is reached, missing the 
correct number by about 0. 1 %. 
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Fig. |26l The approximations appear to converge rapidly. But the limit does not coincide with the known exact value, 
although it happens to lie numerically quite close. Extrapolating the successive approximations by an extremely 
accurate fit to the analytically known large-order behavior [151 with a function ^[fpiateau^'^ ~ 0.843) - A + B L^", we 
find convergence to A = 1 - 0.001 136, which misses the correct Umit A = 1. The other two parameters are fitted best 
by B = -0.002495 and k = 0.922347 (see Fig.|27l). 

We may easily convince ourselves by numerical analysis that the error in the limiting value is indeed linked to 



the failure of the strong-coupling behavior ( 11431 ) to have the power structure of lll5ll . For this purpose we change 



the function f(x) in equation ( 1140b slightly into f{x) — > f(x) = f(x) + 1, which makes the integrals for b,„ in (I141I I 
convergent. The exact limiting value 1 of A remaines unchanged, but b[^^ acquires now the correct strong-coupling 



power structure of II15I1 . For this reason, we can easily verify that the application of variational theory with a dynamical 
determination of co yields the correct strong-coupling limit 1 with the exponentially fast convergence of the successive 
approximations for L — > oo like » 1 - exp(-1.909 - 1.168 L). 

It is worthwhile emphasizing that an escape to complex zeros which BR propose to remedy the problems of the 
5-expansion is really of no help. It has been claimed [53] and repeatedly cited [49], that the study of the anharmonic 
oscillator in quantum mechanics suggests the use of complex extrema to optimize the 5-expansion. In particular, the 
use of so-c ailed /am/Z/e,? of optimal candidates for the variational parameter z has been suggested. We are now going 
to show, that following these suggestions one obtains bad resummation results for the anharmonic oscillator Thus we 
expect such procedures to lead to even worse results in field-theoretic applications. 

In quantum mechanical applications there are no anomalous dimensions in the strong-coupling behavior of the 
energy eigenvalues. The growth parameters a and o) can be directly read off from the Schrodinger equation; they are 
a = 1 /3 and oj = 2/3 for the anharmonic oscillator (see Appendix A). The variational perturbation theory is applicable 
for all couplings strengths g as long as b'^\z) becomes stationary for a certain value of z- For higher orders L it must 
exhibit a well-developed plateau. Within the range of the plateau, various derivatives of b'^\z) with respect to z will 
vanish. In addition there will be complex zeros with small imaginary parts clustering around the plateau. They are, 
however, of limited use for designing an automatized computer program for localizing the position of the plateau. The 
study of several examples shows that plotting b'^\z) for various values of a and oj and judging visually the plateau is 
by far the safest method, showing immediately which values of a and co lead to a well-shaped plateau. 

Let us review briefly the properties of the results obtained from real and complex zeros of d^,b^Q\z) for the anhar- 
monic oscillator In Fig.|28] the logarithmic error of ii^^' is plotted versus the order L. At each order, all zeros of the 
first derivative are exploited. To test the rule suggested in [53], only the real parts of the complex roots have been used 
to evaluate b[^\ The fat points represent the results of real zeros, the thin points stem from the real parts of complex 
zeros. It is readily seen that the real zeros give the better result. Only by chance may a complex zero yield a smaller 
error. Unfortunately, there is no rule to detect these accidental events. Most complex zeros produce large errors. 
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Figure 28: Logarithmic eiTor of the leading strong-couphng coefficient b''^ ^ of the ground state energy of the anharmonic oscillator with potential. 
The errors are plotted over the order L of the variational perturbation expansion. At each order, all zeros of the first derivative have been exploited. 
Only the real parts of the complex roots have been used to evaluate . The fat points show results from real zeros, the smaller points those from 
complex zeros, size is decreasing with distance from real axis. 



We observe the existence of families described in detail in the textbook 161] and rediscovered in Ref. 1531]. These 
families start at about N = 6, 15, 30, 53, respectively. But each family fails to converge to the correct result. Only 
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Figure 29: Deviation of the coefficient from tlie exact value is shown as a function of perturbative order L on a linear scale. As before, fat dots 
represent real zeros. In addition to Fig.|28] the results obtained from zeros of the second deiivative of fojj''* are shown. They give rise to own families 

with smaller eiTors by about 30%. At A' = 6, the upper left plot shows the start of two families belonging to the first and second deiivative of 
respectively. The deviations of both families are negative. On the upper right-hand figure, an enlargement visualizes the next two families starting 
alN = 15. Their deviations are positive. The bottom row shows two more enlargements of families starting at W = 30 and N = 53, respectively. 
The deviations alternate again in sign. 



a sequence of selected members in each family leads to an exponential convergence. Consecutive families alternate 
around the correct result, as can be seen more clearly in a plot of the deviations of b^^^ from their L — > oo -limit in 
Fig. |29] where values derived from the zeros of the second derivative of fej,^' have been included. These give rise to 
accompanying families of similar behavior, deviating with the same sign pattern from the exact result, but lying closer 
to the correct result by about 30%. 



9. Appendix B: Ground-State Energy from Imaginary Part 

We determine the ground state energy function Eoig) for the anharmonic oscillator on the cut, i.e. for g < in the 
bubble region, from the weak coupling coefficients ai of equation (|97] l. The behavior of the ai for large I can be cast 
into the form 

L 
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The Pj can be determined by a high precision fit to the data in the large / region of 250 < I < 300 to be 

-/^ _1 ^ ill 391691 40783 1915121357 10158832895 70884236139235 
' ~2' 24' ~6~' 3456 ' 48 ' 248832 ' 124416 ' 71663616 ' 
60128283463321 286443690892 144343264152266 351954117229 2627843837757582 



4478976 1423 43743 6 2339 

230619387597863 12122186977970425 41831507430222441029 



10 24 3550 

where the rational numbers up to j - 6 are found to be exact, whereas the higher ones are approximations. 



Equation ( 1 1451 1 can be read as recurrence relation for the coefficients a/. Now we construct an ordinary differential 
equation for E{g) :— E'^l^^yig) from this recurrence relation and find: 



E{g) = . (147) 



All coefficients being real, real and imaginary part of E(g) each have to satisfy this equation separately. The point 
g = 0, however, is not a regular point. We are looking for a solution, which is finite when approaching it along the 
negative real axis. Asymptotically E{g) has to satisfy E(g) ^ exp(l/gy6_i) - exp(l/3g). Therefore we solve d 1471 1 
with the ansatz 



E(g) = g" exp 



(148) 



to obtain a - -1/2 and 

95 619 200689 2229541 104587909 7776055955 9339313153349 172713593813181 



'1,2,3,... 



24' 32 ' 1152 ' 1024 ' 3072 ' 12288 ' 688128 ' 524288 

(149) 

1 248602386820060039 1453 1 808399402704 1 603 1 663 1 1257983672027964 1 73696056792 1 



139886592 54391637278720 1435939224158208 

1 0905 1 8247 1 7547897884794645746723 455740 1 7678 1 73074497482074500364087 



348951880031797248 3780312033677803520 



This is in agreement with equation ( llOOI l and an improvement compared to the WKB results of 115211 . Again, the first 
six rational numbers are exact, followed by approximate ones. 

10. Appendix C: First-Order Differential Equations for E„{g) 
Given a one-dimensional quantum system 

(Ho+gV)\n,g)^E„(g)\n,g) (150) 

with Hamiltonian H - Hq + g V , eigenvalues E„(g) and eigenstates \n,g) we consider an infinitesimal increase dg in 
the coupling constant g. The eigenvectors will undergo a small change: 

\n,g + dg)^\n,g) + dg^u„k\k,g) (151) 

SO that 

° k*n 
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Given this, we take the derivative of ( 1150b with respect to g and multiply by (m, ^| from the left to obtain: 



{m, g\V- E'„(g)\n, g)^Yj ^l-^o + ^ ^ " g) ■ (153) 



Setting now m-n and m + n'm turn, we find: 



E'nig) ^Vnnig) (154) 

V,„„(g) =M„,„ {EJg) - E„(g)) , (155) 

where Vm„(g) = {m,g\V\n,g). 

Equation (I154l l governs the behavior of the eigenvalues as functions of the coupling constant g. In order to have a 
complete system of differential equations, we must also determine how the V„,n{g) change, when g changes. With the 
help of equations (1152b and ( 1155b . we obtain: 



v: 



VmkVkn V VmkVkn 



Z ymkykn ymkykn nsi\ 

ki^m ki^n 

Equations ( 11541 ) and ( 1157b together describe a complete set of differential equations for the energy eigenvalues E„{g) 
and the matrix-elements V„m{g). The latter determine via ( 11551 ) the expansion coefficients u„,„(g). Initial conditions 
are given by the eigenvalues E„(Q) and the matrix elements V„m{0) of the unperturbed system. 

[I] J. Goldstone, A. Salam, and S. Weinberg, Phys. Rev. 127, 965 (1962). 
[2] C. De Dominicis. J. Math. Phys. 3, 983 (1962). 

[3] R.P. Feynman, Statistical Mechanics, Benjamin, Reading, 1972. 

[4] R.P. Feynman and H. Kleinert, Phys. Rev. A 34, 1986 (http//physik. fu-berlin.de/-klelnert/159/159.pdf). 
[5] H. Kleinert Phys. Lett. A 173, 332 (1993) (http://physik.fu-berlin.de/~kleinert/213 i. 

[6] H. Kleinert, Path Integrals in Quantum Mechanics. Statistics, Polymer Physics, and Financial Markets, World Scientific, Singapore, 2009 

(http : //www.physlk. f u-berlin. de/"'klelnert/b5 I. 
[7] H. Kleinert, Annals of Physics 266, 135 (1998) (http://physik.fu-berlln.de/~kleinert/255/255.pdf I. 

[8] J.A. Lipa, D.R. Swanson, J.A. Nissen, T.C.P Chui and U.E. Israelsson, Phys. Rev. Lett. — bf 76, 944 (1996); J.A. Lipa, D. R. Swanson, J.A. 

Nissen, Z.K. Geng, PR. Williamson, D.A. Strieker, T.C.P Chui, U.E. Israelsson and M. Larson, Phys. Rev. Lett. — bf 84, 4894 (2000). 
[9] H. Kleinert, Phys. Rev. D 60, 085001 (1999) (hep-th/9812197 I (hep-th/9812197 1; Phys.Lett. A 277, 205 (2000) (cond-mat/9906107 I. 
[10] PM. Stevenson, Phys. Rev. D 30, 1712 (1985); D 32, 1389 (1985); PM. Stevenson and R. Tarrach, Phys. Lett. B 176, 436 (1986). 

[II] W. Janke, H. Kleinert, Phys. Rev. Lett. 75, 2787 (1995) (quant-ph/9502019 1. 

[12] E.J. Weniger, Phys. Rev. Lett. 77, 2859 (1996); see also F Vinette and J. Ci'zek, J. Math. Phys. 32, 3392 (1991). 
[13] H Kleinert and W. Janke, Phys. Lett. A 206, 283 (1995) ( quant-ph/9509005 I. 

[14] R. Guida, K. Konishi, H. Suzuki, Ann. Phys. (N.Y.) 249, 109 (1996). 

[15] H. Kleinert, Phys. Rev. D 57, 2264 (1998) (cond-mat/9801167 I; Phys. Rev. D 58, 107702 (1998) (co nd-mal/9803268t . 
[16] B. Hamprecht and H. Kleinert, Phys. Rev. D 68, 065001 (2003) ( hep-th/03021 16 1. 

[17] R. Seznec and J. Zinn-Justin, J. Math. Phys. 20, 1398 (1979). This paper has developed important techniques for understanding the conver- 
gence mechanism of vaiiational perturbation theory. 

[18] H. Kleinert and V. Schulte-Frohlinde, Critical Phenomena in O'^-Theory, World Scientific, Singapore, 2001 
(htt p : / /www.physlk . f u-berlin. de/"kleinert/b8 I. 

[19] J. Zinn-Justin, Quantum Field Theory and Critical Phenomena, Clarendon, Oxford, 1989. 

[20] See Section 19.4 in [H]. 

[21] RJ. Wegner, Phys. Rev. B 5, 4529 (1972); B 6, 1891 (1972). 

[22] B.G. Nickel, D.L Meiron, and G.B. Baker, Univ. of Guelph preprint 1977 (unpublished). The preprint is readable on the WWW at 
|http: / /www.physlk . f u-berlin. de/"'kleinert/nickel/guelph.pdf 
The results are cited and used in Chapters 1 9 and 20. They were axtended to seven loops by 
D.B. Murray and B.G. Nickel, Univ. of Guelph preprint 1991. 

The additional coefficients of the renormalization group functions are listed in Section 20.4 of the textbook fl^ . 

[23] See Section 20.2 in 1 18]. 

[24] H. Kleinert, Phys. Rev. D 60, 085001 (1999) ftep-th/9812197». See Fig. 9. 

[25] M. Holzmann, G. Baym, J.-P Blaizot and F Laloe, Phys. Rev. Lett. 87, 120403 (2001). 

[26] R Arnold, G. Moore and B. Tomasik, cond-mat/0107124 

[27] J. Zinn-Justin, Quantum Field Theory and Critical Phenomena (Oxford University Press, Oxford, England, 1996). 

[28] M.BijlsmaandH.T.C.Stoof, Phys. Rev A 54, 5085 (1996) 

[29] FF de Souza Cruz, M.B. Pinto and R.O. Ramos, Phys. Rev. A 65, 053613 (2002) ( cond-mal/0112306) . 



31 



[30] E. Braaten, and E. Radescu, (cond-mat/0206186|. 

[31] These are the results of f^cf]. They diflfer from those of Refs. f3] which are 03 =0.644519, «4i =0.87339, 042 = 3.15905, 043 = 1.70959, 044 = 
4.4411,045 = 2.37741. The coefficients of the series )57t to be resummed differ mainly in the last term: /_i = -126.651 10"*, /o = 0,/i = 
-4.04857 10""*, /2 = 2.40587 10"'', /3 = -2.06849 10"''. 

[32] H. Kleinert, Strong-Coupling Behavior of Phi'^ -Theories and Critical Exponents, Phys. Rev. D 57 , 2264 (1998); Addendum: Phys. Rev. D 
58 , 107702 (1998) (cond-mat/9803268 I; Seven Loop Critical Exponents from Strong-Coupling (/>'*-Theory in Three Dimensions, Phys. Rev. D 
60 , 085001 (1999) (hep-th/98 12197 1; Theory and Satellite Experiment on Critical Exponent alpha of Specific Heat in Superfiuid Helium 
Phys. Lett. A 277, 205 (2000) ( cond-mat/9906107 I. 

[33] H. Kleinert, Strong-Coupling (p'^-Theory in 4 - e Dimensions, and Cri tical Exponent, Phy s. Lett. B 434 , 74 (1998) tcond-mat/9801167t ; 
Critical Exponents without beta-Function, Phys. Lett. B 463, 69 (1999) l icond-maV'9906359) . 

[34] See (Hl^lil and references cited there. 

[35] With standard normalization conditions used in the 3-dimensional (/■''-theory, the approach to scaling is governed by Wegner's exponent a) 
(see l33l ). The present definition of m differs from the inverse correlation length m = f "' by a factor: m = mZ^' oc m m"'''- for m — ► 0. This 
changes the exponent of approach to cj' = w/(l - rj/l). I thank B. Kastening for noting this. 

[36] P Arnold and G. Moore, Phys. Rev. Lett. 87, 120401 (2001); Phys. Rev. E 64, 066113 (2001). The authors derive a l/N correction factor 
(1 - 0.521 IN) to the leading iV ^ 00 result. 

[37] V.A. Kashurnikov, N.V. Prokof'ev and B.V. Svistunov, Phys. Rev. Lett. 87, 120402 (2001). 

[38] G. Baym, J.-P Blaizot M. Holzmann, F. Laloe and D. Vautherin, Phys. Rev. Lett. 83, 1703 (1999). 

[39] G. Baym, J.-P Blaizot and J. Zinn-Justin, Europhys. Lett. 49, 150 (2000). 

[40] P. Arnold and B. Tomasik, Phys. Rev. A62, 063604 (2000). This paper starts out from the 3+1-dimensional initial theory and deiives from it 
the three-dimensional eft'ective classical field theory, the field-theoretic generalization of the quantum-mechanical eft'ective classical potential 
of R.P Feynman and H. Kleinert, Phys. Rev. A 34, 5080 (1986). This reduction program was started for the Bose-Einstein gas by A.M.J. 
Schakel, Int. J. Mod. Phys. B 8, 2021 (1994); J. Mod. Phys. B 8, 2021 (1994); Boulevard of Broken Symmetries, Habihtationsschrift, FU- 
Berlin, I cond-mat/9805 1 52 1 (1998). Unfortunately, Schakel did not go beyond the one-loop level so that he was happy to have found a positive 
shift AT^/Tc, and did see the cancellation at the two-loop level. See his recent paper in J. Phys. Stud. 7, 140 (2003) (cond-mat/030105Cl) . 

[41] FF de Souza Cruz, M.B. Pinto and R.O. Ramos, Phys. Rev. B 64, 014515 (2001). 

[42] P Grueter, D. Ceperley, F Laloe, Phys. Rev. Lett. 79, 3549 (1997) ( cond-mat/9707028 I. 

[43] See Eq. (20.23) in the textbook [18] or S.E. Derkachov, J.A. Gracey, and A.N. Manashov, Eur. Phys. J. C 2, 569 (1998) ftep-ph/9705268». 

[44] W. Janke and H. Kleinert, Phys. Lett. A 117, 353 (1986) http : //www .physik. f u-berlin. de/~kleinert/133 1; Phys. Rev. Lett. 58, 144 
(1986). H. Kleinert, Phys. Lett. A 257, 269 (1999) I cond-mat/98 1 1 308 1; M. Bachmann, H. Kleinert, A. Pelster, Phys. Lett. A 261, 127 (1999) 
lcond-mat/9905397 1; Physical Review E 63, 051709/1-10 (2001) lcond-mat/001 1281 1; see also B. Kastening, Phys.Rev. A 68, 061601 (2003) 
I cond-mat/0303486 1; Phys.Rev. A 69, 043613 (2004) I cond-mal/03090 60 1; Phys. Rev. E 73, 01 1 101 (2006) (cond-mal/0508614) ; Phys. Rev. 
A 70, 043621 (2004) ( cond-mat/0406035 1. ' ' 

[45] J.S. Langer, Ann. Phys. 41, 108 (1967). 

[46] CM. Bender and T.T. Wu, Phys. Rev. 184, 1231 (1969); 

[47] S. Coleman, Nucl. Phys. B 298, 178 (1988). 

[48] B.Bellet, PGarcia, A.Neveu, Int. J. of Mod. Phys. A 11, 5587(1997) 
[49] J.-L. Kneur, D. Reynaud, lhep-th/0205133v2). See also (sOl, (Hi, Q. 
[50] E. Braaten, E. Radescu, lcond-math/0206186 1'l). 
[51] The low-order results were first obtained by 

H. Kleinert, Phys. Lett. B 300, 261 (1993) Ihttp : / /www . physik ■ f u-berlln ■ de/-kleinert/214) , 

and extended by 

R. Karrlein and H. Kleinert, Phys. Lett. A 187, 133 (1994) lhep-th/9504048». 
[52] J. Zinn-Justin, J. Math Phys. 22(3), 511 (1981). The first 10 coefiicients of expansion (99) are calculated. 

[53] B.Bellet, P.Garcia, A.Neveu, Int. J. of Mod. Phys. All, 5587(1997). The family structure of optimal variational parameters emphasized in 
this paper was discussed in great detail eariier in Chapter 5 of the textbook [6], but with correct application rules. 

[54] M. Bachmann, H. Kleinert, and A. Pelster, Phys. Rev. A 62, 52509 (2000) lquant-ph/0005074 1, Phys. Lett. A 279, 23 (2001) (quant- 
ph/000510). 

[55] J.E. Avron, B.G. Adams, J. Cizek, M. Clay, M.L. Glasser, P Otto, J. Paldus, and E. Vrscay, Phys. Rev. Lett. 43, 691 (1979). 
[56] L.D. Landau and E.M. Lifschitz, Quantenmechanik, Sechste Auflage (Akademie-Verlag BerHn, 1979). 

[57] J.E. Avron, I.W. Herbst, B. Simon, Phys. Rev. A 20, 2287 (1979). See also the resummation treatment in J.-C. Le Guillou, J. Zinn-Justin, 
Ann. Phys. 147, 57 (1983). 

[58] H. Ruder, G. Wunner, H. Herold, and F. Geyer, Atoms in Strong Magnetic Fields (Springer- Verlag, Berlin, 1994). 
[59] G. Ahlers, Phys. Rev. A 3, 696 (1971); K.H. Mueller, G. Ahlers, F Pobell, Phys. Rev. B 14, 2096 (1976); 

[60] L.S. Goldner, N. Mulders and G. Ahlers, J. Low Temp.Phys. 93 (1992) 131. 

[61] R. Guida and J. Zinn-Justin, Critical exponents of the N-vector model, J. Phys. A 31, 8130 (1998) (cond-mal/9803240 > 

[62] J.C. Le Guillou and J. Zinn-Justin, Phys. Rev. Lett. 39, 95 (1977); Phys. Rev. B 21, 3976 (1980); J. de Phys. Lett 46, L137 (1985). 

[63] H. Kleinert and V. Schulte-Frohlinde, J. Phys. A 34, 1037 (2001) (cond-mat/9907214 1. 

[64] D.B. Murray and B.G. Nickel, unpubHshed. 

[65] A. Pellisetto and E. Vicari, preprint IFUP-T H 52/97,|cond-mat/9711078| 
[66] F. Jasch and H. Kleinert, BerHn preprint 199 (cond-inat/9907214) . 
[67] W. Janke, Phys.Lett. A148 (1990) 306. 

[68] H.G. Ballesteros, L.A. Fernandez, V. Martin-Mayor and A. Munoz Sudupe, Phys. Lett. B 387, 125 (1996). 
[69] M. Ferer, M.A. Moore, and M. Wortis, Phys. Rev. B 65, 2668 (1972). 
[70] P Butera and M. Comi, Phys. Rev. B 56, 8212 (1997) (hep-Iat/9703018t . 



32 



